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The recoil correction and spin-orbit force for the possible B*B* and D*D* states 
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In the framework of the one-boson exchange model, we have calculated the effective potentials between two 
heavy mesons B*B* and D*D* from the t- and u-channel tt-, 77 -, p-, a>- and cr-meson exchanges. We keep the 
recoil corrections to the and D*5* systems up to O(^), which turns out to be important for the very loosely 
bound molecular states. Our numerical results show that the momentum-related corrections are favorable to the 
formation of the molecular states in the = L, 7^^ = in the B*B* and D'D* systems. 
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I. INTRODUCTION 


A lot of charmonium-like states have been reported in the 
past decade by the experiment collaborations such as Belle, 
Barbar, CDF, DO, LHCb, BESIII, and CLEOc. The under¬ 
lying structures of many charmonium-like states are not very 
clear. Sometimes they are called as XYZ states. They de¬ 
cay into conventional chromium, but not all of them can be 
accommodated into the quark-model charmonium spectrum. 
The neutral XYZ states include X(3872) H, T(426m H, 
T(4008) il, T(4360) lH, T(4660) lH, and T(4630) B| etc. 
There are also many charged charmonium-like states such as 
Zi(4050) andZ2(4250) 1^^(4485) HU, 2^(3900) CMl, 
Zc(4020) lfl3l] . Zc(4025)1[m 1- The charged bottomonium-like 
states Zi(10610) and Zi(10650) were observed by Belle Col¬ 
laboration B. 

Theoretical speculations of these XYZ states include the 
hybrid meson ifl^ . tetraquark states dynamically 

generated resonance ll24l] and molecular states etc. 

Since many of these XYZ states are close to the thresholds 
of a pair of charmed or bottom mesons, the molecular hypoth¬ 
esis seems a natural picture for some of these states. 

Within the framework of the molecular states, there exist 
extensive investigations of the charged Z^ and Z* states & 
S. In our previous work JUt], we explored the possibility 
of Zc(3900) as the isovector molecule partner of X(3872) and 
considered the recoil correction and the spin-orbit force in the 
DD* and BB* system very carefully. 

Although there exist quite a few literatures on the possibil¬ 
ity of Zc(4025) as the D*D* molecular state and Zi(10650) as 
the B^B* molecular state, most of the available investigations 
are either based on the heavy quark spin-flavor symmetry or 
derived in the mg —» 00 . In other words, the recoil correc¬ 
tion and the spin-orbit interaction have not been investigated 
for the D*D* and B*B* systems. Since the binding energies of 
these system are very small, the high order recoil correction 
and the spin-orbit interaction may lead to significant correc¬ 
tions. 
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In this work, we will go one step further. We will consider 
the recoil correction and the spin-orbit force for the D*D* and 
B*B* systems. With the one-boson-exchange model (OBE), 
we will derive the effective potential with the relativistic La- 
grangian and keep the momentum related terms explicitly in 
order to derive the recoil correction and the spin-orbit interac¬ 
tion up to (9(1/M“), where M is the mass of the heavy meson. 
We investigate the B*B* system with = \^,J^ - l‘'“ for 
Zfc(10650), and the D*D* system with - l‘'“ 

for Zc(4025). For completeness, we also investigate the B*B* 
system and D*D* system with other quantum numbers: — 

\-J^ = 1 ++, = 1 ++, and Q = I'-. 

Compared to the DD* case, the expressions of the recoil cor¬ 
rections and spin orbit force are more complicated. There ap¬ 
pear several new structures. For some systems, the numerical 
results show that the high order correction is important for the 
loosely bound heavy-meson states. 

This paper is organized as follows. We first introduce the 
formalism of the derivation of the effective potential in Sec. 
In] We present our numerical results in Sec. HII] and Sec. m 
The last section is the summary and discussion. 


II. THE EFFECTIVE POTENTIAL 


A. Wave function, Effective Lagrangian and Coupling 
constants 


First, we construct the flavor wave functions of the isovec¬ 
tor and isoscalar molecular states composed of the B*B* and 
D*D* as in Refs, [j^ . The flavor wave function of the 
B*B* system reads 

' |1,1) = |B*+B*''), 

. |1,-1) = |B*^B*°), (1) 

|1,0) = ■j=(\B*^B*-) - \B*°B*°)), 


|0,0) = —(\B*^B*-) + |B*'’b*°» (2) 

V2 

For the D*D* system 

' |1,1) = |D*°D*+), 

, |1,-1) = |D*0d*-), (3) 

|1,0) = ■^(\D*°D*°) - |D*-D*+», 
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|0,0) = 4=(|D*°D*°) + \D*-D*^)) (4) 

V2 

The meson exchange Feynman diagrams for the B*B* and 
D*D* systems at the tree level is shown in Fig. [T] 


B'.D" B*,D* 

-*-^ - * - 

I 

I 

g* ,5* _ \ _ B\ £)• 


FIG. 1: The Feynman diagrams for both the D'D* and 6*B* systems 
at the tree level. 


Based on the chiral symmetry, the Lagrangian for the pseu¬ 
doscalar, scalar and vector meson interaction with the heavy 
flavor mesons reads 


■Cp - 


-il±MPl^d^cPbaP\ + i^MPi,d^c/>baPf 

Jn Jn 

/tt Jn 


Zp = -i^-MP%<PabPT-i^MPfd^,(P,hPh 


fn 


fn 


+ ^dPpfd‘^4,,bP*^e,,ap - yPfd‘^‘f>ahd>^Pl''€pyafii6) 

Jn Jn 


-Cv - 


i^-^PKMl-i^d,PbVlPl 

- iJ2Ag„e^apy&^Pbd‘^vlr:^ 

- iJ2Ag,e,,f,,Pl^d‘‘vld''Pl 

q p*t I r)*V]/P p*t 

' ^^b ^ba^t^^va ’■ ^ba^va 

- il j2Ag,M*Pl^{d^Vy - dyVf.haPl'^, 


(7) 


£v = -i^d,Plv:^Pb + i^Plvy,Pb 

+ iJ2Ag„e^apyP7^d^vl^d''Pb 
+ i'^-^d,^XbJ 7 -i^^Xbdpn^ 

- i2 j2Ag,M*PT\d^Vy - dyV.Ulf, 


( 8 ) 


£s = -2g,MPb(TPl + 2g,M*P;VF;; (9) 


Jls = -2gMPl(rPa + 2g,M*P*JaCTPT ( 10 ) 

where the heavy flavor meson fields P and P* represent P = 
(D°,D+) or (B-,B°) and P* = (D*°,.D*+) 01 ^( 5 *-, . 8 *°). Jts 
corresponding heavy anti-meson fields P and P* represent P = 
{D°,D-) or ( 8 +, 8 °) and P* = (D*°,D*-) or ( 8 *+, 8 *‘’). (j), V 
represent the the exchanged pseudoscalar and vector meson 
matrices, cr is the only scalar meson interacting with the heavy 
flavor meson. 


(p = 


71° , n 

V 2 V 6 


71° , Jl 

V 2 V 6 1 


V = 


P° I OJ 

V 2 V 2 
P" 


P° I OJ 

^/2 Vl ‘ 


( 11 ) 


( 12 ) 


According to the OBE model, five mesons ( tt, cr, p, oj and 
77 ) contribute to the effective potential. For the D*D* and B*B* 
systems, the potentials are the same for the three isovector 
states in Eqs. (1)~(4) with the exact isospin symmetry. Ex¬ 
panding the Lagrangian densities in Eqs. (5)~(10) leads to 
each meson’s contribution for this channel. These channel- 
dependent coefficients are listed in Table |T] 

The pionic coupling constant g - 0.59 is extracted from the 
width of /,r=132 MeV is the pion decay constant. 

According the vector meson dominance mechanism, the pa¬ 
rameters gy and p can be determined as g,, = 5.8 and p - 0.9. 
At the same time, by matching the form factor obtained from 
the light cone sum rule and that calculated from the lattice 
QCD, we can get A - 0.56 GeV ' J^S^. The coupling 
constant related to the scalar meson exchange is gs — g„l2 V6 
with g;r = 3.73 1^ . All these parameters are listed in 

TableHIl 


TABLE I: coefficients 



isospin 

meson-exchange 

P 

OJ 

cr 

n 

'7 

ZTO* 

/ = 1 

-1/2 

1/2 

1 

- 1/2 

1/6 

/ = 0 

3/2 

1/2 

1 

3/2 

1/6 

B*B* 

/ = 1 

-1/2 

1/2 

1 

- 1/2 

1/6 

/ = 0 

3/2 

1/2 

1 

3/2 

1/6 


In order to include all the momentum-related terms in our 
calculation, we introduce the polarization vector of the vector 
mesons. At the rest frame we have 


Q = ( 0 , el) 


(13) 


We make a lorentz boost to Eq. [13] to derive the polarization 
vector in the laboratory frame 


lab 


= (- 


P(P- el) 
m(PQ + m) 


m 


(14) 
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TABLE II: The coupling constants and masses of the heavy mesons 
and the exchanged light mesons used in our calculation. The masses 
of the mesons are taken from the PDG 1^ 



mass(MeV) 

coupling constants 

pseudoscalar 

= 134.98 
m, = 547.85 

g = 0.59 

A = 132MeV 

vector 

mp = 775.49 
= 782.65 

gv = 5.8 
/? = 0.9 

A = 0.56GeV-' 

scalar 

77V = 600 

gj =g,r/2V6 

gn = 3.73 

heavy flavor 

miy, - 2010.25 
mB> - 5325.0 



where p - {po, p) is the particle’s 4-momentum in the labora¬ 
tory frame and m is the mass of the particle. 

B. Effective potential 

With the wave function and Feynman diagram, we can de¬ 
rive the relativistic scattering amplitude at the tree level 


instead of p' and p in the practical calculation. 

In the OBE model, a form factor is introduced at each vertex 
to suppress the high momentum contribution. We take the 
conventional form for the form factor as in the Bonn potential 
model. 


^ A^-ml ^ A^-ml 
A 2-^2 A 2 + ^ 


( 21 ) 


nia is the mass of the exchanged meson and m* is the mass of 
the heavy flavor meson D* or B*. So far, the effective potential 
is derived in the momentum space. In order to solve the time 
independent Schrodinger equation in the coordinate space, we 
need to make the Fourier transformation to V{q, k). The details 
of the Fourier transformations are presented in the Appendix. 

The expressions of the potential through exchanging the cr, 
p mesons are 


Va- = -Ca-gliib ■ eJXet' ■ 

~^o-g~s ^[(^^ ' ^<3 )(^6' ' )F3rio- + S 

+Co-g]-^i{eb' ■ x et) ■ L]F 5 ,o- (22) 

m ^ 


</|5 ID = + i{f\T\i) = 6 fi + (2nf5\pf - pdiMp, (15) 

where the T-matrix is the interaction part of the S-matrix and 
Mfi is defined as the invariant matrix element. After apply¬ 
ing Bonn approximation to the Lippmann-Schwinger equa¬ 
tion, the S-matrix reads 


{f\S\i) = 5fi - 2n6(Ef - EdiVp (16) 

with Vfi being the effective potential. Considering the differ¬ 
ent normalization conventions used for the scattering ampli¬ 
tude Mfi, T -matrix Tfi and Vfi, we have 


02 2 

-Cp^iet ■ eV)(et' ■ ej ')Fi,p 

-Cp2A^gl(4 X ej)(et' x ej ^)F 2 tp 
+Cp2A^gl[(eb X ej)(eb> x ej ')F 3 ,ip -H S i2F3t2p] 
Bhl 2A/3gl 


-Cpi 


2 ot *2 

q2„2 


m 


-)[(ft • ea)(eb' ■ ej )F3,ip + S i2F3,2p] 


1 . 


-Cp^(4 ■ e«^)(e?' ■ e:/)[F4np + {--V^F4,2p)] 

■ eJ'')[(f«^' X et) ■ L]F5,p (23) 

2 ot *2 m* 


Vfi--' 


M 


fi 


n 2 / 7 /° n 2 p,' 


M 


fi 


n ' 2 .m /° n 2 ,mP 


(17) 


where ppf denotes the four momentum of the final (initial) 
state. 

During our calculation. Pi = (E\,p) and P 2 = (E 2 ,—p) 
denote the four momenta of the initial states in the center mass 
system, while P 3 = (£ 3 , 77 ') and P 4 = (£ 4 ,-/ 7 ') denote the 
four momenta of the final states, respectively. 

^ = P3-Pl =(£3 -£i,p"'-/?) = (£2-£4,^ (18) 


The oj and p meson exchange potentials have the same form 
except that the meson mass and channel-dependent coeffi¬ 
cients are different. 

The expression of the potential through exchanging the n 
meson is 

2 

= C,^[(etxe^^)(e?-xej')£3,i. + 5i2£3,2j (24) 
Jtt 

where S j 2 and S 12 have the form 

^'12 = [3(F- et)(r ■ e^^) - (4 ■ ■ ej^) (25) 


is the transferred four momentum or the four momentum of 
the meson propagator. For convenience, we always use 

q^p'-p (19) 

and 

^=2^ip’+P) ( 20 ) 


5 12 = 3[r ■ (4 X efl''')][r X (et' ■ ej^)] - (4, x 4j)(eb' x ej ^026) 

Compared to the DD* case, there appear several new in¬ 
teraction operators: (4b x 4a)(4' x Sa’^), 8 ^ 2 , S 12 and i( 4 ’ ■ 
6 /^)[(^^ X ^) ■ Z^]. These operator represent the new form of 
the tensor, spin-spin and spin-orbit interactions. 
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Similarly, the i] and n meson exchange potential has the 
same form in the D*D* and B*B* system except the meson 
mass and channel-dependent coefficients. The explicit forms 
of Tf,ta,T^,ua,T^tva, T^uva are shown in the Appendix. 

In our calculation, we explicitly consider the external mo¬ 
mentum of the initial and final states. Due to the recoil cor¬ 
rections, several new terms appear which were omitted in the 
heavy quark symmetry limit. These momentum dependent 
terms are related to the momentum k - j{p' + p): 


P 

(p- +m^ 


(27) 


and 


/(ej ■ ea’^)[(ej x et) -(kxq)] 


-f) . 9 

q- + mi 


(28) 


The term in Eq. ( l28l l is the well-known spin orbit force. In 
short, all the terms in the effective potentials in the form of 
Fstip, F 4 iip, F^tp etc with the sub-indices 3,4,5 arise from the 
recoil corrections and vanish when the heavy meson mass m* 
goes to infinity. The recoil correction and the spin orbit force 
appear at (9(1/M^). 


C. Schrodinger equation 


With the effective potential V(F) in Eqs. (23) ~ (27), we are 
able to study the binding property of the system by solving the 
Schrodinger Equation 

(_|!v2 + y(i^-£)'I'(i^ = 0, (29) 

2p 

where 'I'(i^ is the total wave function of the system. The total 
spin of the system 5=1 and the orbital angular momenta 
L - 0 and L - 2. Thus the wave function 'l'(i^ should have 
the following form 


'P(^ = <As(^ + M4 (30) 


where (As(fO and are the S -wave and D-wave functions, 
respectively. We use the same matrix method in Ref. ll45l] to 
solve this S-D wave couple-channel equation. 

We detach the terms related to the kinetic-energy-operator 
from V(F) and re-write Eq. ( l29b as 

(-^v2-^[VV0 + a(0V'] 

2p 2p 

-i-y(i^ - £ )'T(f) = 0 (31) 


with 



in which a{r) is 


(32) 


a{r) 


o 2 2 

{-2p)[-Cp^{P, ■ eV){et ■ ei,^)rAap 

a2 2 

-Cp|||(et • eVXet' ■ eJ^)r4,2^] (33) 


The total Hamiltonian contains three angular momentum 
related operators (et • eJXet' ■ e?'^), (el, x eJXet' x ej"^), 
S 12 , i(Fh' ■ ea’^)[iPi^ x ^b) ■ L], which corresponds to the spin- 
spin interaction, spin orbit force and tensor force respectively. 
They act on the S and D-wave coupled wave functions +4>d 
and split the total effective potential V(f) into the subpotentials 
yss(F), Vsnir), Vosif) and VoDir)- The matrix form reads 


{4>s + (t>D\ (4 ■ ePXet' ■ ei4)V(F) \4>s + 4 >d) 
_ I Vss(f) 0 ) 

\ 0 Vooif) ) 


(34) 


{4>s + 4 >d\S \ 2 ^(r)\(ps + ^d) - 


y-^yosir) 


_\ 

1 


(35) 


{(ps + (f>D\ (4 X ePXet’ X ei-^)V(r) |0s + (po) 
^(Vss(r) 0 \ 

\ 0 Vooir)) 


{<Ps + MS X 2 y{p\<ps +4 >d)-[_ ^yM 


{4>s + M ii4’ ■ efl'^)[(efl' X 4,) ■ pyiP \^s + ^d) 
_/0 0 \ 
io IVDD(r)) 


(36) 



(38) 


III. NUMERICAL RESULTS EOR THE 5*B* SYSTEM 


We diagonalize the Hamiltonian matrix to obtain the eigen¬ 
value and eigenvector. If there exists a negative eigenvalue, 
there exists a bound state. The corresponding eigenvector is 
the wave function. We use the variation principle to solve the 
equation. We change the variable parameter to get the lowest 
eigenvalue. We also change the number of the basis functions 
to reach a stable result. 


A. Z,,( 10650) 

Since the mass of the charged bottomonium-like state 
Zi,( 10650) is close to the B*B* system, we first consider the 
possibility of the B*B* molecule with 7*^ = 1+, - l'*'“. In 

order to reflect the recoil correction of the momentum-related 
terms, we plot the effective potential of the S-wave and D- 
wave with or without the momentum-related terms in Eig. |2] 
Vs and Vd are the effective potentials of the S -wave and D- 
wave interactions after adding the momentum-related terms, 
y' and y^ are the effective potentials of the S -wave and D- 
wave interactions without the momentum-related terms. Eig. 
|2]C corresponds to the - 1+, - l+“ B*B* system, 

where the curves of y^ and y', and Vd and V'^ are almost over¬ 
lapping. In other words, the recoil correction is small. 
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We collect the numerical results in Table |III] E and E' 
are the eigenenergy of Hamiltonian with and without the 
momentum-related terms, respectively. The fourth, fifth and 
sixth column represent the contribution of S -wave, D-wave, 
and spin-orbit force components, respectively. The last col¬ 
umn is the mass of 8 *8* as a molecular state of = 1^, 
- l"^". When the cut off lies within 2.2 - 2.8 GeV, there 
exists a bound-state. The binding energy with the recoil cor¬ 
rection is between 0.97-15.15 MeV. The binding energy with¬ 
out the recoil correction is between 0.94 - 14.98 MeV. When 
the cutoff parameter A = 2.2 GeV, the binding energy is 0.97 
MeV, and the recoil correction is only -0.03 MeV. The con¬ 
tribution from the spin-orbit force is as small as 0.001 MeV. 
When the cutoff parameter A = 2.8 GeV, the binding energy is 
15.15 MeV, and the recoil correction is -0.07 MeV. The corre¬ 
spondence spin-orbit force contribution is 0.02 MeV, which is 
also small compared with the binding energy. The recoil cor¬ 
rection and the contribution of the spin-orbit are very small. 
However the recoil correction is favorable to the formation of 
the 8 *8* molecular state with 7*^ = 1+, - l'^". 

TABLE III: The bound state solution of the 8*8* system with 7° = 
1 ^(in unit of MeV) and different A. E and £' is the eigen¬ 
energy of the system with and without the momentum-related terms 
respectively. We also list the separate contribution to the energy from 
the S-wave, D-wave and spin-orbit force components respectively in 
the fourth, fifth and sixth column. The last column is the mass of the 
8*8* system with 7*^ = H, J''‘~ = as a molecular state. 


A(GeV) 



Eigenvalue 


Mass 


total 

S 

D 

LS 

(MeV) 

2.2 

E 

-0.97 

13.15 

0.15 

0.001 

10649.03 

E' 

-0.94 

-12.91 

-0.15 

- 

10649.06 

2.4 

E 

-3.49 

-28.68 

0.32 

0.004 

10646.51 

E' 

-3.43 

-28.27 

0.32 

- 

10646.57 

2.6 

E 

-8.04 

-49.81 

0.56 

0.01 

10641.96 

E' 

-7.94 

-49.19 

0.55 

- 

10642.06 

2.8 

E 

-15.15 

-77.53 

0.88 

0.02 

10634.85 

E' 

-14.98 

-76.66 

0.87 

- 

10635.02 


From Fig[3]C, it is clear that the n exchange is much more 
important than the other meson-exchanges. Considering that 
the coupling constant g is extracted from the D* decay width 
with some uncertainty, we multiply g by a factor from 0.99 
to 1.1 to check the dependence of the binding energy on this 
parameter. The numerical results are listed in Table HV] The 
binding energy with the recoil correction varies from 6.39 - 
36.81 MeV. The binding energy without the recoil correction 
varies from 6.3 - 36.57 MeV. The binding energy is sensitive 
to the coupling constant. 

B. The 8 * 8 * system with 7° = C, 

For the 8 *8* system with 7*^ = 1, = 1^^, there is no 

bound state when the cutoff varies in a reasonable range. We 
multiply the coupling constant g by a factor to investigate the 


TABLE IV: The 8*8* system with 7® = H, J'‘‘' = (in units of 
MeV) and different coupling constant g and A = 2.0 GeV. The other 
notations are the same as in Tablelllll 


A(GeV) 



Eigenvalue 


Mass 


Total 

5 

D 

LS 

(MeV) 

g ■ 0.99 

E 

-6.39 

-43.28 

0.49 

0.01 

10643.61 

E' 

-6.30 

-42.70 

0.49 

- 

10643.7 

8 

E 

-8.04 

-49.81 

0.56 

0.01 

10641.96 

E' 

-7.94 

-49.19 

0.55 

- 

10642.06 

g-l.Ol 

E 

-9.91 

-56.68 

0.63 

0.01 

10640.09 

E' 

-9.79 

-56.02 

0.63 

- 

10640.21 

g-l.l 

E 

-36.81 

-132.98 

1.32 

0.02 

10613.19 

E’ 

-36.57 

-131.96 

1.31 

- 

10613.43 


dependence of g. We list the results in Table IVl The cutoff pa¬ 
rameter is fixed at A = 2.0 GeV. There appears a bound state 
when g increases by a factor 1.3. When the factor changes 
from 1.3 - 1.5, the binding energy with recoil correction is 
between 0.82 - 19.25 MeV. The binding energy without the 
recoil correction is between 0.8 - 19.18 MeV. When the bind¬ 
ing energy is 0.82, the recoil correction is -0.02 MeV, and the 
contribution of the spin-orbit is 0.004 MeV. When the binding 
energy is 19.25 MeV, the the recoil correction is -0.07 MeV, 
and the contribution of the spin-orbit force is 0.04 MeV. From 
Fig|2]D, the effect of the recoil correction is very small for the 
8 *8* system with - 1, - 1^+. But it is favorable 

for the formation of the molecular state. Fig[3]D shows the 
contributions of each meson-exchange. 


TABLE V: The 8*8* system with 7® = ff, J‘‘‘- = (in units of 
MeV) with the variation of the coupling constant g and A = 2.0 GeV. 
The other notations are the same as in Tablelllll 


A(GeV) 



Eigenvalue 


Mass 


Total 

5 

D 

LS 

(MeV) 

g-L3 

E 

-0.82 

1.99 

-2.23 

0.004 

10649.18 

E' 

-0.80 

1.97 

-2.21 

- 

10649.2 

g-L4 

E 

-6.68 

l.Ti 

-8.94 

0.02 

10643.32 

E' 

-6.63 

7.69 

-8.92 

- 

10643.37 

g-L5 

E 

-19.25 

15.41 

-18.71 

0.04 

10630.75 

E' 

-19.18 

15.38 

-18.70 

- 

10630.82 


C. The 8*8* system with = 0 , J‘‘‘^ = H 

We also investigate the 8 *8* system with = 0“, = 

l^V The effective potential and the meson-exchange contri¬ 
bution are shown in Fig and [3] A. The recoil correction is 
large. The numerical results are listed in Table I VII The cutoff 
varies from 1.0- 1.3 GeV. The binding energy with the recoil 
correction is around 0 - 41.69 MeV. The binding energy with¬ 
out the recoil correction is between 0.3 - 53.47 MeV. When 
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FIG. 2: The effective potential of the S*B* system. The Labels 
A,B,C, and D correspond to the four cases I = 0“, = 1^“; / = 0^, 

jpc ^ ^++. j ^ jpc ^ J+-. j ^ jpc ^ j++^ respectively 
from top to bottom. and Vd are the effective potentials of the 
S -wave and D-wave interactions with the momentum-related terms, 
while V' and Vj are the S -wave and £l-wave effective potentials with¬ 
out the momentum-related terms. 


FIG. 3: The effective potential from the different meson exchange in 
the 5*B* system. Labels A,B,C,D are the same as in Fig. |2] 


D-wave contribution. But it is unfavorable for the formation 
of the molecule. 


the cutoff parameter A = 1.0 GeV, there is no bound state af¬ 
ter adding the recoil correction. The recoil correction is 0.3 
MeV. When the cutoff parameter A = 1.3 GeV, the binding 
energy is 41.69 MeV. The recoil correction is 12.05 MeV, and 
the contribution of the spin-orbit force is 4.3 MeV. In this case, 
the recoil correction is significant and it is almost as big as the 


D. The S*B* system with = 0^, 

For the B*B* system with = 0+, = 1++, V, and 

are clearly different within the range of 0 - 0.6 fm as can be 
seen from Fig[2]B. Fig[3]B shows the contributions of each 
meson-exchange. The recoil correction in this system is im- 
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FIG. 4: The effective potential of the D*£)* system. The Labels 
A,B,C, and D correspond to the four cases I = 0“, = 1^“; / = 0^, 

jpc ^ ^++. j ^ j + ^ jpc ^ J+-. j ^ jpc ^ j++^ respectively 
from top to bottom. and Vd are the effective potentials of the 
S -wave and D-wplws interactions with the momentum-related terms, 
while V' and Vj are the S -wave and £l-wave effective potentials with¬ 
out the momentum-related terms. 


portant. The results are listed in Table IVIII When the vari¬ 
ation of the cutoff parameter is between 0.9 -1.2 GeV, the 
binding energy with the recoil correction changes from 6.71 
to 59.74 MeV. The binding energy without the recoil correc¬ 
tion changes from 6.86 - 63.76 MeV. When the cutoff param¬ 
eter A = 0.9 GeV, the binding energy is 6.71 MeV. The re¬ 


FIG. 5: The effective potential from the different meson exchange in 
the 0*5* system. Labels A,B,C,D are the same as in Fig.|4] 


coil correction is 0.15 MeV, and the contribution of the spin- 
orbit force is 0.03 MeV. When the cutoff parameter A = 1.2 
GeV, the binding energy is 63.76 MeV. The recoil correction is 
4.02 MeV, and the contribution of the spin-orbit force is 1.39 
MeV. The recoil correction is of the same order as the D-wave 
contribution. Moreover, it increase with the binding energy. 
Therefore, the recoil effect can not be neglected. However, 
the recoil correction is unfavorable for the formation of the 
molecule in this system. 
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TABLE VI: The bound state solutions of the B'B* system with 7® = 
0“, J‘‘'~ = (in unit of MeV) with the cutoff A. The other notations 
are the same as in Tablelllll 


A(GeV) 



Eigenvalue 


Mass 


total 

S 

D 

LS 

(MeV) 

1.0 

E 

- 


- 


10650 

E' 

-0.30 

-0.59 

-0.07 

- 

10649.7 

1.1 

E 

-4.34 

-1.88 

-3.08 

0.42 

10645.66 

E' 

-5.51 

-2.19 

-4.29 

- 

10644.49 

1.2 

E 

-17.07 

-9.05 

-7.92 

1.68 

10632.93 

E' 

-21.79 

-12.09 

-11.12 

- 

10628.21 

1.3 

E 

-41.69 

-24.61 

-14.80 

4.30 

10608.31 

E' 

-53.74 

-30.81 

-24.77 

- 

10596.26 


TABLE VII: The bound state solutions of the B'B" system with 7° = 
0^, = 1^^ (in unit of MeV) with the cutoff A. The other notations 

are the same as in Tablelllll 


A(GeV) 



Eigenvalue 


Mass 


total 

S 

D 

LS 

(MeV) 

0.9 

E 

-6.71 

-22.27 

0.18 

0.03 

10643.29 

E' 

-6.86 

-22.85 

0.15 

- 

10643.14 

1.0 

E 

-19.36 

-43.77 

0.25 

0.18 

10630.64 

E' 

-20.15 

-46.10 

0.03 

- 

10629.85 

1.1 

E 

-37.51 

-65.91 

0.53 

0.58 

10612.49 

E' 

-39.58 

-71.13 

-0.2 

- 

10610.42 

1.2 

E 

-59.74 

-86.12 

1.34 

1.39 

10590.26 

E' 

-63J6 

-94.99 

-0.43 

- 

10586.24 


IV. NUMERICAL RESULTS EOR 77*5* SYSTEM 
A. 77*5* system with 7® = U, J’’'- = 

Due to the isospin symmetry, the interaction in the D*D* 
system has the same form with the B*B* system. There¬ 
fore, we repeat the same investigations for the D*D* system. 
Zc(4025) was observed in the recoil mass spectrum in the 
process — > {D*D*)^7t^ ifTil] . The mass of Zc(4025) is 

close to the threshold of D*D*, and Zc(4025) has the quantum 
with - 1^, - l"^". Therefore we first consider the 

possibility of the D*D* system as the molecular state with the 
quantum number 7^^ = 1^, = l^L 

From Fig|4]C, the curves of Vj and V' are almost overlap¬ 
ping, same as those of Vd and V^, which indicates that the re¬ 
coil correction is small. Fig|5]C shows that n meson-exchange 
plays an important role in the interaction. 

Unfortunately, we did not get a bound state within a reason¬ 
able range of the cutoff parameter and coupling constant. The 
value of the pionic coupling constant was extracted from the 
decay width of the D* meson where the pion is on the mass- 
shell. However we need the value of the coupling constant in 
the potential where the pion is off-shell. Considering the big 


influence of the n meson-exchange, we multiply the coupling 
constant g by a factor to check the dependence of the results 
on g. The cutoff parameter is fixed at A = 2.0 GeV. The results 
are shown in Table IVIlTl 

When the factor reaches 1.6, there appears the bound state. 
The binding energy with the recoil correction is 1.15 MeV 
and the recoil correction is -0.05 MeV. The contribution of 
the spin-orbit force is 0.01 MeV. When the factor is 1.8, the 
binding energy with recoil correction is 28.62 MeV, and the 
recoil correction is -0.28 MeV. The contribution of the spin- 
orbit force is 0.07 MeV. As in the B*B* system with - H, 
= l^A the recoil correction is not so big. But the recoil 
correction is favorable for the formation of D*D* molecular 
state. In other words, the existence of the D*D* molecule de¬ 
pends on the coupling constant g sensitively. 


TABLE VIII: The 77*£»* system with 7^ = 1+, = i+- (in units 

of MeV) with the variation of coupling constant g and A = 2.0 GeV. 
The other notations are the same as in Tablelllll 


A(GeV) 



Eigenvalue 


Mass 


T otal 

S 

77 

LS 

(MeV) 

g- 1.6 

E 

-1.15 

-22.96 

0.37 

0.01 

4019.35 

E' 

-1.10 

-22.27 

0.36 

- 

4019.40 

g-1.7 

E 

-10.19 

-75.21 

1.10 

0.04 

4010.31 

E’ 

-10.03 

-74.02 

1.08 

- 

4010.47 

g-1.8 

E 

-28.62 

-139.80 

1.80 

0.07 

3991.88 

E' 

-28.34 

-138.16 

1.77 

- 

3992.16 


B. 77*D* system with 7® = 1 , 

We also investigate the D*D* system with - 1“, = 

1^^. From Fig|4]D and Figl^D, the recoil correction is very 
small while the n meson-exchange plays an significant role in 
the interaction. There also does not exist a bound state when 
the cutoff parameter is within a reasonable range. We also 
study the variation with the coupling constant g. The results 
are shown in Table IIXI 

The recoil correction is also very small in this system. For 
example, when the binding energy is 4.71 MeV, the recoil cor¬ 
rection is only -0.03 MeV, and the contribution of the spin- 
orbit force is 0.05 MeV. And the recoil correction is favorable 
for the formation of the molecular state. 


C. The 77*D* system with 7° = 0 , 7^^ = U and 7® = 0*^, 

JPC ^ 1 ++ 

For the isoscalar D*D* system, and V' are very different 
from Fig. |4]A and B. Fig|5]A and B show the contributions 
of each meson-exchange. The recoil contribution is large and 
unfavorable for the formation of the molecular states. There 
exist bound states for the isoscalar D*D* system with different 
C-parity. The results are listed in Tables 11^ and IXII 
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TABLE IX: The £)*D* system with 7® = 1, 7^'- = 1++ (in units of 
MeV) with the variation of coupling constant g and A = 2.0 GeV. 
The other notations are the same as in Tablelllll 


A(GeV) 



Eigenvalue 


Mass 


T otal 

S 

D 

LS 

(MeV) 

^•2.4 

E 

-4.71 

7.72 

-8.47 

0.05 

4015.79 

E' 

-4.68 

7.69 

-8.49 

- 

4015.82 

^■2.5 

E 

-11.05 

13.22 

-14.96 

0.09 

4009.45 

E' 

-11.02 

13.19 

-15.02 

- 

4009.48 

^•2.6 

E 

-20.37 

19.46 

722.76 

0.13 

4000.13 

E' 

-20.34 

19.43 

-22.86 

- 

4000.16 


For the case, when cutoff parameter changes 

from 1.4 - 1.6 GeV, the binding energy with the recoil correc¬ 
tion is within 0.58 - 17.09 MeV. The binding energy without 
the recoil correction is within 3.83 - 17.09 MeV. The recoil 
correction is large. For example, when the binding energy is 
0.58 MeV, the recoil correction is 3.25 MeV, even bigger than 
the binding energy itself. The contribution of the spin-orbit 
force is 0.4 MeV, which is almost as large as the binding en¬ 
ergy. 

For the = 1^^ case, when cutoff parameter changes 
from 1.3 - 1.6 GeV, the binding energy with the recoil correc¬ 
tion is within 9.13-43.25 MeV. The binding energy without 
the recoil correction is within 10.59 - 49.23 MeV. The recoil 
correction is also significant. For example, when the bind¬ 
ing energy is 9.13 MeV, the recoil correction is 1.46 MeV. 
The contribution of the spin-orbit force is 0.75 MeV, which 
is almost as large as D-wave contribution. When the binding 
energy is 43.25 MeV, the recoil correction is 5.98 MeV. The 
contribution of the spin-orbit force is 5.21 MeV. 


TABLE X: The D'D* system with 7° = 0 , J'’‘' = 1^ (in unit of 
MeV). The other notations are the same as in Tablelllll 


HGeV) 



Eigenvalue 


Mass 


total 

S 

D 

LS 

(MeV) 

1.4 

E 

-0.58 

-5.03 

-0.43 

0.40 

4019.92 

E’ 

-3.83 

-15.14 

-3.06 

- 

4016.67 

1.5 

E 

-5.67 

-19.23 

-1.45 

1.91 

4014.83 

E' 

-17.25 

-44.41 

-8.84 

- 

4003.25 

1.6 

E 

-17.09 

-40.12 

-2.52 

4.76 

4003.41 

E' 

-42.61 

-79.61 

-17.73 

- 

3977.89 


V. SUMMARY AND DISCUSSION 

With the one-boson-exchange model, we have systemati¬ 
cally studied the possible loosely bound B*B* and D*D* sys¬ 
tems with (1) = O'", = 1++, (2) = 0-, = 1+-, 

(3) = 1+, = r- and (4) = r, = 1++. We 

consider the n,ria-,p and oj meson exchange in the derivation 


TABLE XI: The 77*7)* system with 7** = 0*^, (in unit of 

MeV). The other notations are the same as in Tablelllll 


A(GeV) 



Eigenvalue 


Mass 


total 

S 

D 

LS 

(MeV) 

1.3 

E 

-9.13 

-34.43 

1.04 

0.75 

4011.37 

E' 

-10.59 

-42.22 

-0.15 

- 

4009.91 

1.4 

E 

-18.20 

-49.42 

2.25 

1.66 

4002.3 

E' 

-20.87 

-61.74 

-0.29 

- 

3999.63 

1.5 

E 

-29.64 

-63.30 

4.22 

3.11 

3990.86 

E' 

-33.82 

-80.59 

0.57 

- 

3986.68 

1.6 

E 

-43.25 

-76.05 

7.17 

5.21 

3977.25 

E' 

-49.23 

-98.50 

1.10 

- 

3971.27 


of the potential. We keep the momentum dependent terms in 
the polarization vector of the two heavy mesons and introduce 
the momentum-related terms in the interaction, which lead to 
the recoil correction and spin-orbit force at 0(1 IM^). 

The B*B* system with - 1+, = l"^" can form the 

molecular state, which may correspond to the heavier Z/, state 
observed by Belle collaboration. When the cutoff parameter 
is within 2.2 - 2.6 GeV, the binding energy is between 0.97 - 
15.15MeV. The recoil correction is small. The contribution of 
the spin-orbit force is also very small. For example, when the 
binding energy is 15.15 MeV, the recoil correction is -0.17 
MeV. The contribution of the spin-orbit force is 0.02 MeV. 
But the recoil correction is favorable to the formation of the 
molecular state. On the other hand, our results shows that the 
binding energy is sensitive to the pionic coupling constant. 

For the isoscalar B*B* system, there exist a bound state 
when changing the cutoff parameter. For the 7^** = 1^“ state, 
when cutoff parameter is within 1.0 - 1.3 GeV, the binding 
energy is between 0-41.69 MeV. For the 7^** = 1^^ state, 
when cutoff parameter is within 0.9 -1.2 GeV, the binding 
energy is between 6.71 - 59.74 MeV. The recoil correction of 
the two systems are both large and important. However, they 
are unfavorable to the formation of the molecular states. 

For the D*D* system with - H, 7^** = l"^", we are 
unable to obtain the bound state within a reasonable cutoff 
range and the pionic coupling constant g extracted from the 
D* decay width. If we enlarge the pionic coupling constant 
by a factor of 1.6 - 1.8, there appears the bound state with the 
binding energy around 1.15-28.62 MeV. The recoil correction 
is small. For example, when the binding energy is 28.62 MeV, 
the recoil correction is -0.28 MeV. The contribution of the 
spin-orbit force is 0.07 MeV. The recoil correction is favorable 
for the formation of the molecular state. 

For the isoscalar D*D* system, there exist bound states 
when changing the cutoff parameter. For the state, 

when the cutoff parameter is within 1.4 - 1.6 GeV, the binding 
energy is around 0.58 - 42.61 MeV. For the 7^*^ = 1^^ state, 
when the cutoff parameter is within 1.3 - 1.6 GeV, the binding 
energy is around 9.13 -43.25 MeV. The recoil correction is 
significant but unfavorable to the formation of the molecular 
states. 
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VI. APPENDIX 


We collect the lengthy formulae in the appendix. 

. exp(mc,r) 

Y(mar) = —;- 


(39) 


ZOhar) = (1 + + -—^—^)Y{mar) (40) 

m^r {niarY 


J 3 ta — it .9 , 1 " 


+ cp- 


+ ml 


= -cn ■ cf2[ mlAY(Ar) - mlYim^r) 

+ (A^ - ml)A^^ ] 

+ IrS\2[-mlZ{mar) + A^Z{Ar) 


Ziifhar) - {-— + —— :^)Y(mar) (41) 

m^r (mary 


+ (A^ - ml)(l + Ar)^F(Ar)] 
= (cTi • cf* i)'F2t\ + S n'Fiti 


(51) 


sin(m„r) 3 sin(m„r) 1 cos(mar) 

Z {mar) = -;--;-+ —-^-;-. (42) 

mar mar mar {mar)^ mar 

1 sin(m„r) 1 cosifhar) 

Zpmar) = -;-+ ——^;- 

niar mar {mar)'^ mar 

where for system DD* 

- ml. 


= fnl,p,a,,n - ('«D - 


for system BB* 




- (mg - ms) . 


(43) 

(44) 

(45) 

(46) 


J2>ua — U 7 o ^ 


A^ + p- 


-O , - 9 


= 3*^' ■ (^2\fnlAY{Ar) - mlYiniar) 
+ {A^-ml)A^—] 

+ ^5 n[-ma^Z{Tnar) + A^Z(Ar) 

+ (A2-m^)(l+Ar)|F(Ar)] 

= (CTI • (T l)J^2u\a + S l2'F2u2a 


(52) 


T\ta - 


ISp — ml 1 
T{{-^^)- 


A'^ + p p + ml 


= m„y(m„r) - AF(Ar) - (A - m„)^^ (47) 


. A^ - ml ^ {p ■ pip ■ p , 
r J- = I 


^lua — 


A^ — ml 1 


A^ + p p + ini' 


-Ar 


- fiiaYiniar) - AY{Ar) - (A - ( 4 ^) 


T2ta = T{{ 


A^ - ml p 


A^ + p p + ml 


- ml[AY{Ar) - maY{mar)] 

+ (A^-mPA^-Y- 


T 2 ua = T{{ 


A^ - ml p 


A^ +p p + ml 


- ml[AY{Ar) - niaYipr)] 


-Ar 


+ {A^-mpA— 


(49) 


(50) 


+ (A^-m^)A—] 

+ ^5i2[m^Z'(mar) +A^Z(Ar) 

+ {A^-ml)(l+Ar)^Y{Ar)] 
- (P ' (TjJT'sula + S \2'F'2u2a 


TAta = m 
A 


A^ 


A^ + p p + ml 


tfi~' 

= ^F(m,r)-—y(Ar) 


A^ - ml . Ar , e 


-Ar 


i-T-P - 

^2 r 

- ^{V^, niaYiniar) - AF(Ar) ■ 

1 9 

- t^4/lQ- + {-2^ P4l2a} 


(53) 


A^-ml 


A 


(54) 
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TAua - Tii 




A2 + ^ ^ + ml 


fh^ A^ 

= -^Yithar)- —Y{Ar) 


A^ - ml , Ar 


e 


-Ar 


2 r 


- \maY{mar)-AY{Ar)- 

- T^ii\a + {--;F4«2q-) 


A^ — m; 


2 ^-Ar 


A 


(55) 


rzrr _ ^ (^1 ' ^)(^2 ' , 


^ 

A^ + + ml 

- ^ ■ z![-OT^Zi(m„r) -H A^Zi(Ar) 

+ [A^-mlf-^] 


- ^ ■ LTstOo 


^5ua — ~7 ^ 

A^ + ^ f- + ml 
- S ■ L[-mlZi{mar) + A^Zi(Ar) 


+ (A-" 




— ^ ■ IZT'^uda 


(56) 


- (A2-m2)A^] 

1/c ^ 3 COS(OT<,r) 2 

+ -(5 12 + cFi • Cr2)[ (-2- -H -:-;-)'«„ 

i mar mar mar 

- (1 + |-)(A)2F(Ar) - (A2 - ml)CA + -)^ ]V 

Ar r 2A 

1 ^ ^ cos(OT„r) 

- ^(S 12 + tr-’i ■ (f2)[ma : - A^Ar) 

3 mo-r 

„-kr 


rT] ■ 0~2 ^ / o ^ ^ \rT~t 

— ^ A 6wlo' ^ ' ^A)A 6 h2(2' 

1 ^ ^ , 

- T^iS\2 + O-1 ■ O-2)T 6uia 


(60) 


(57) 


r,„. = 


A? + ^ (p- + i 


= 5'-L[- 


mf,Zj(mar) -I- A^Zi(Ar) 
-Ar 


+ (A^-mi)—] 
- S ■ LT'suOa 


(58) 


^ -mlpcri-k)(cZ2-k), 

y 6ua — -^iVrT -tn y 


r 


[P- +ml 


-[m„y(m„r)-(A)^y(Ar) 


(A-_,„2)a_] 


+ :r(5i2 + CTi • cr2)[ (1 + 2 :—)OT„y(Ar) 

3 mar 

- (1 + i-)(A)2y(Ar) 

Ar 

7 p-Ar 

- {A^-ml)(A+-)^]V 

- ^( 5 12 + CTi ■ (f'2)[maY{mar) - AY(Ar) 

- (A= -»4)i^iv= 

cr [ ■ ( 7^2 1 ^ ^ 

~ ^ P^6u\a "t" 12+^1 ■ 0'2)P^6u2a 

- 2 (‘^ 12 + cfi ■ A l)T6u3a 
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